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1 Introduction 

The magnetoencephalogram (MEG) observes the 
brain activity as a magnetic field by the intracellular 
current. The time resolution of the method is better 
than 1 ms. However, it is difficult to estimate 
multiple sources and spread sources, because the 
Equivalent Current Dipole method (ECD) is 
generally used for the source analysis of MEG. 

On the other hand, the functional Magnetic 
Resonance Imaging (fMRI) observes the changing 
of the cerebral blood flow of local area. Though it is 
inferior in time resolution, it is possible to observe 
spread sources. Therefore, MEG and fMRI have 
complementary roles in the brain research [1]. 
Recently, Weighting Minimum Norm Estimation 
method (w-MNE) based on Minimum Norm 
Estimation method (MNE) has been proposed 
[2][3]. In this paper, we incorporate fMRI data into 
weighting of w-MNE. First we measure MEG and 
fMRI on same subject for same experimental tasks. 
Second the proposed analysis method for the MEG 
and fMRI data is examined. 

2 Method 

The measurement of fMRI and estimation of the 
location of the activity area in the brain are 
performed in prior to the measurement of MEG. 

Let us define the magnetic field measured by the /th 
sensor as b i? and a vector B = [bi b 2 * * -b m ] T is 
defined as a set of the measured data. Here, m is the 
total number of the sensors and superscript T 
denotes the matrix transpose. Uniformly distributed 
lattice of n points are assumed in the brain and a 
vector P = [pi P2 ' ' 'Pn] T is defined so that a 
component of source moment for each lattice points. 
If the positions of the source are defined, the 
relation between B and P can be written as the 
linear system of the form: 

B = LP (1) 

where L is an mxn lead field matrix. Generally L 
can be evaluated by the equation of spherical model. 


Our purpose on the biomagnetic inverse problem is 
to calculate P from B and L: 

P = L*B (2) 

where L* is called the Moore and Penrose 
generalized inverse and is an n*m matrix, 

L*=(L T L)' 1 L T . (3) 

In this problem, the generalized inverse matrix of L 
cannot be used, because the solution of L* includes 
noise in B. To avoid this instability, we have used 
the regularized inverse matrix of Wc, which is well 
known as the Thikhonov regulation. Then the 
solution is written as 

P = (L t L + yWc)' 1 L T B (4) 

where Wc and y, respectively, an nxn weighing 
matrix and a scalar regularization parameter. Then 
Wc satisfies 

Wc = diag(|l 1 |, |1 2 |, • • *,|l n |) (5) 

where f is the /th column vector of the lead field 
matrix L. 

In this report, we use a weighting matrix W mri , 
which is incorporating the measurement of fMRI, as 
Wc. Matrix W mri is diagonal, whose elements are 
fundamentally equal to the corresponding elements 
of W c , except for the elements that are being active 
in the fMRI measurement. The value of such 
exceptional elements are forced to be smaller than 
that of the corresponding elements of W c . Then the 
solution is written as 

P = (L t L+ W mri ) ‘L t B . (6) 

Our aim is to minimize the cost function: 

Ec = (B - LP) t (B - LP) + y P T W mri T W mri P (7) 

where the first term means the difference of the 
magnetic fields between the measured data and the 
calculated data, and the second term means the total 
strength of the weighted moment. 



We use a regularization parameter y to adjust the 
balance between the first term and the second term. 
The suitable value of y is determined by the 
computer simulations. 

3 The measurement condition 

3.1 Measurement of MEG 

The MEG system has 65 detection sites for vector 
components. Each detection site contains three 
gradiometers, and three gradiometers are 
perpendicular to each other. Thus, the system 
consists of 195 dc-SQUID first-order gradiometers 
(baseline 50 mm). The typical system noise in the 
magnetic shielded room is 10 fT/VHz. 

The stimulation is the right visual field flashing 
stimulation of 2 Hz. The condition of the data 
acquisition is as follows. (1) Sampling frequency: 
400 Hz, (2) Analog band-pass filter : 0.1-100 Hz, 
notch filter: 50 Hz, (3) Digital band-pass filter 
0.5-40 Hz, Baseline -100-0 ms, (4) Measurement 
time: 6 minutes, Averaging: 230 times. 

3.2 Measurement of fMRI 

The fMRI measurements were carried out with a 
1.5T whole body scanner (STRATIS2: Hitachi 
Medical Co. Japan.), using a standard imaging head 
coil. The subject head was fixed with a head holder 
to minimize movements. The fMRI was performed 
using blood oxygen level dependent (BOLD) 
contrast with a 2D FLASH sequence (echo time TE 
= 45ms, repetition time TR = 4600ms). Six fMRI 
slices around the V1 area in the occipital lobe were 
recorded simultaneously (slice thickness 6 mm, gap 
1mm, voxel size 4.0x4.0x6.0 mm, sagittal slice ). 
During fMRI measurement, resting (12 images) and 
active state (12 images) alternated every 72 s 
(6 s per image). For rest and stimulation, 6 cycles 
each were performed. The visual stimulation were 
presented on the right visual field flashing 
stimulation of 2 Hz. The activation threshold of P < 
0.005 for single voxels was chosen. 

4 The measuring results 

4.1 MEG resulta 

The measurement result of Visual Evoked Field 
(VEF) by the right visual field flashing stimulation 
is shown in Figs. 1-3. Fig. 1 shows the VEF 
waveforms. Figures 2 and 3 show the isomagnetic 
contour map at latency 92 ms and at latency 130 ms, 
respectively. The viewpoint of Figs. 2 and 3 is the 
back of the head. Figs. 2(a) and 3(a) shows normal 


component and Figs. 2(b) and 3(b) shows tangential 
component. Figs. 2 and 3 both show that there was a 
source that seemed to be in the left occipital area 
near the middle. However the estimated location of 
the signal source was in the right occipital area near 
the middle at latency 130 ms. The estimated source 
by ECD is shown in Fig.7(a-l),(b-l). 
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Figure 1: VEF waveforms by the right visual field 


flashing stimulation. 
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Figure 2: Visual evoked field (VEF) observed from 
back of head. The data was at about 92 ms after the 
stimulus by the right visual field flashing 
stimulation, (a) Normal component (b) Tangential 
component. The solid line shows outer component, 
broken lines shows inner component in (a). 
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Figure 3: VEF at about 130 ms after the stimulus, 
(a) Normal component, (b) Tangential component. 
Other conditions are the same as Fig. 2. 
















4.2 fMRI results 

The measuring result of fMRI by the right visual 
field flashing stimulation is shown in Fig. 4. Figure 
4 shows around visual area by 6 slices. The right 
side of Fig. 4 is the left side for the subject, and 
then, the left side is the right for the subject. There 
was a spot that seems to be in the VI area in Fig. 
4(a). The activities near the scalp were the noise. 



Figure 4: The estimated activations of fMRI by the 
right visual field flashing stimulation. The 
viewpoints of figures are from back of the head. The 
activation threshold of P < 0.005 for single voxels 
was chosen. The result of using W mri is shown in 
section 5. 
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5 The application to data of VEF by the 
right visual field flashing stimulation 

The experiment of fMRI has been done and the 
activity area was estimated before MEG experiment. 
W^i was weighted in 10 lattice points. The number 
of lattice points were about 3000 in the whole head 
(lattice point interval: lcm). The MEG data at 92 ms 
and 130 ms were analyzed. Fig. 5 and 6 show the 
arrow map of the sources estimated by Wc and Wmri. 
From Fig. 5(b), estimated signal source by Wmri is 
more local than that estimated by Wc. There are no 
weighting which used prior information of fMRI on 
the point in Fig. 6(b). However the estimated current 
distribution is similar to that estimated by Wc. 

Fig. 7 shows the estimated results using ECD and 
Wmri superimposed on the MR imaging. Fig. 7(a- 
2),(b-2) show the point where the maximum volume 
of the norm for the estimated source in all points. 
Table 1 shows estimated location of the source 
using ECD and Wmri, and it is shown by the head 
coordinate that is determined by outer cranium 
point. From Fig. 7, the position of the signal source 
estimated by ECD and Wmri came out in the V1 area 
at 92 ms. 


Table 1 

'.The signal source using ECD and W mri . 

Latency 

ECD 

Wmi-MNE 

[ms] 

[cm] 

[cm] 

92 

(-1.00 -6.86 5.83) 

(-1 -5 8) 

130 

(1.27 -6.65 7.73) 
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Figure 5: The distribution of the source using W c and W mri for the VEF at 92 ms. (a) The result was estimated 
by W c . (b) The result was estimated by W mri . 
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Figure 5: The distribution of the source using W c c and W mri for the VEF at 130 ms. (a) The result was 
estimated by W c . (b) The result was estimated by W mri . 



































6 Discussion 

We proposed that the technique of signal source 
estimation method using the prior information of f- 
MRI, and it apply to the analysis of VEF to the right 
visual field flashing stimulation. 

In Fig. 7, there is the estimated source position by 
respectively using ECD and Wmn in the V1 area, but 
the position is slightly different. 

The signal source of ECD is in the right-brain at 
130ms. However, in the experimental result of f- 
MRI in Fig. 4, there are no activities in the right- 
brain. In such a case that a contradiction with the 
results of MEG and f-MRI occurs. The estimated 
current distribution by We and Wm are almost the 
same (see Fig. 6). The fact means that the estimation 
by does not cause any bad effect even in such a 
contradiction occurs. 

The premise of this analysis is "Estimated result of 
MEG and estimated result of f-MRI are almost 
same”. Further the remedy except for it must be 
examined. In future we will consider not only the 
weighting using the prior information on locations 
of f-MRI and but also weighting using the blood 
mobilization rate. 
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Figure 7: The result of estimated by ECD at (a-1) 92ms and (b-l) 130ms, respectively. The 
result of estimated by W mri at (a-2) 92ms and (b-2) 130ms, respectively. 









